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Abstract 

We compute decay constants of heavy-light mesons in quenched lattice QCD with a lattice 
spacing of a ~ 0.04 fm using non-perturbatively 0(a) improved Wilson fermions and 
0(a) improved currents. We obtain f Ds = 220(6)(5)(11) MeV, f D = 206(6)(3)(22) MeV, 
f Bs = 205(7) (26) (17) MeV and f B = 190(8) (23) (25) MeV, using the Sommer parameter 
ro = 0.5 fm to set the scale. The first error is statistical, the second systematic and 
the third from assuming a ±10% uncertainty in the experimental value of ro- A detailed 
discussion is given in the text. We also present results for the meson decay constants fx 
and f n and the p meson mass. 



PACS: 12.38.Gc, 13.20.Fc, 13.20.He 



1. Weak decays of heavy-light mesons with c and b quarks are interesting for studies of CP 
violation and determination of the CKM mixing angles. New experimental data on such decays 
are emerging (e.g. [1-4]) and their interpretation requires knowledge of hadronic matrix elements 
governed by the strong interaction. Lattice QCD allows one to calculate the strong matrix 
elements from first principles. However, if the heavy quark mass tuq is of the order of the 
inverse lattice spacing a, considerable discretization effects proportional to powers of arriQ 
occur. 

One possibility for coping with this problem is to use an effective theory such as Heavy Quark 
Effective Theory (HQET) [5] or Nonrelativistic QCD (NRQCD) [6]. These formalisms start 
from an infinitely heavy quark and consider corrections to this limit in the form of an expansion 
in the inverse of tuq. However, to study the charm quark in HQET or NRQCD requires a 
considerable number of correction terms, and one still has to worry about the uncertainty from 
the truncation of the 1 / tuq expansion. A formulation for relativistic quarks where masses can 
be of 0(1) in lattice units is the Fermilab approach [7] and modifications thereof, developed 
in [8,9] and [10]. One expects that the dominating discretization effects are then proportional 
to powers of momenta of O(Aqcd)- While within HQET non-perturbative renormalization is 
possible [11], in many of the calculations using effective theories the renormalization constants 
are calculated only in perturbation theory (e.g. in Ref. [12]), leading to further uncertainties. 

Another possibility is to simulate on very fine lattices, and this is the approach we have 
adopted in the present paper. We have performed a quenched lattice study of heavy mesons 
with a lattice spacing a of about 0.04 fm. On such a fine lattice a relativistic treatment of the 
charm quark should be justified and we expect that discretization errors are small compared to 
previous calculations on coarser lattices. We also make an attempt to study B mesons in our 
relativistic framework. Even on our fine lattice we cannot simulate B mesons directly, but the 
required extrapolation becomes relatively short-range. We expect that the resulting uncertainty 
is not much larger than the systematic error caused by the use of an effective theory. For 
example, recent unquenched calculations of fs and fs s [12,13] employ NRQCD for the b quark 
and quote a ~ 10% error based on perturbation theory and other systematic effects. 

In this article we present results for the leptonic decay constants of the D s , B s , D and B 
mesons. We also evaluate light meson masses and decay constants to compare with previous 
quenched calculations of the light spectrum on coarser lattices and in order to be able to dis- 
entangle discretization and quenching effects. 

2. Our results are based on the analysis of 114 quenched Wilson gauge configurations simulated 
at the coupling parameter (3 = 6.6 with a mixed heatbath and microcanonical overrelaxation 
algorithm using the publicly available MILC code [14]. The lattice volume is 40 3 x 80, i.e. our 
lattice extends over 40 points (~ 1.59 fm) in space and 80 points in time. The lattice spacing 
is determined using the Sommer parameter ro = 0.5 fm. This choice is motivated by a previous 
calculation [15] which used r to determine the lattice spacings and found that the results for 
fo s from a quenched lattice and a lattice with Nf = 2 agreed (a rs 0.1 fm in these calculations). 
From the interpolating formula given in [16], one finds for our lattice a" 1 = 4.97 GeV. 

For the quarks we use the 0(a) improved clover formulation [17], with the nonperturbative 
value of the clover coefficient cgw — 1-467 determined in Ref. [18]. We work with seven quark 
masses corresponding to three "light" hopping parameters k = 0.13519,0.13498,0.13472 and 
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four "heavy" hopping parameters, k = 0.13000,0.12900,0.12100,0.11500. Statistical errors are 
estimated by means of a bootstrap procedure using 500 bootstrap samples. For the central 
values we take the median. The error bars are calculated including 34% of the sample values 
below and above the median, respectively. Since the upper and lower error bars are found to be 
quite symmetric for most of our data, we just quote the larger of the two. The autocorrelation 
times for the pseudoscalar meson propagator appear to be small. In the worst case we studied, 
the autocorrelations decay after a distance of one configuration. 

To extract the decay constants we follow the procedure described in Ref. [19]. For light and 
for heavy-light mesons we calculate the correlation functions 

X 

c s Mt) = vj2(n^t)P s \o)), (i) 

X 

where A 4 is the local axial vector current operator, P the pseudoscalar density which can be 
local (i = L) or Jacobi smeared (i = S), and V is the spatial lattice volume. 
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0.13519 


0.13519 


0.1059(13) 


0.1928(62) 


0.0376(11) 


0.0248(13) 


0.0312(09) 


0.13498 


0.13519 


0.1231(12) 


0.2005(48) 


0.0388(11) 


0.0251(12) 


0.0324(09) 


0.13498 


0.13498 


0.1388(10) 


0.2107(43) 


0.0403(10) 


0.0258(10) 


0.0337(08) 


0.13472 


0.13519 


0.1422(11) 


0.2065(39) 


0.0405(11) 


0.0261(11) 


0.0339(09) 


0.13472 


0.13498 


0.1560(10) 


0.2186(33) 


0.0418(10) 


0.0270(10) 


0.0351(08) 


0.13472 


0.13472 


0.1722(09) 


0.2292(27) 


0.0434(10) 


0.0283(09) 


0.0366(08) 



Table 1: Light meson masses and decay constants in lattice units. 
Masses and amplitudes are determined from fits of the correlation functions with 

C s PP (t) = Af P (e- Et + e- E ^) , (2) 
C S AS) = A s P L A4 (e- Et -e- E ^), (3) 

where E is the ground state energy. In Table [TJ we give the raw data for the light pseudoscalar 
meson masses, determined from Cpp, and for light vector meson masses from smeared-local 
correlation functions of the spatial components of the vector currents. 

To determine the bare quark masses, we calculate K cr u, the k value corresponding to massless 
quarks, from a fit of the squared mass of a pseudoscalar meson ( "pion" ) consisting of quarks 
with mass parameters Ki and k 2 as a function of the averaged 0(a) improved quark mass 

(amps) 2 = aiarhg, (4) 

with 

rh q = (1 + b m am q ) m q , m q = -{m ql + m q2 ) 

and am q i = |(^- — ), i = 1,2. We use the non-perturbative value of —0.6636 for the 
improvement parameter b m using an interpolating formula from Ref. [20]. The fit includes all 
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data with > 0.13472, where we find the improved quark masses to just lie on a straight 
line. We find K cr u = 0.135472(11). The hopping parameter corresponding to the average u and 
d quark mass, Kg, is determined by setting nips on the left hand side of Eq. (jl]) equal to the 
physical pion mass, mps = 138 MeV. We find K£ = 0.135456(10). 

We parameterize the quark mass dependence of light meson decay matrix elements with 
hopping parameters K\ and K2 by fitting them to a function of the form 

c + ci afh q . (5) 

The light quark mass dependence of masses and decay matrix elements of heavy-light mesons 
is parameterized using a linear fit as in Eq. (jSJ), with rh q being the light quark mass instead of 
the average quark mass. 

We also calculate the vector ( "p meson" ) mass. The fit and the chiral extrapolation assuming 
a quark mass dependence as in Eq. (jHJ) are shown in Fig.fTJ At ni we find 846(37) MeV (the error 
is statistical), which is roughly a 10% (2a) discrepancy with experiment. We compare our result 
to other recent quenched calculations in Table [2j Within errors our result agrees with Ref. [21], 
where a continuum extrapolation from coarser lattices with 0(a) improved clover fermions is 
performed. We also list studies employing chiral lattice fermions where smaller quark masses 
can be reached while coarser lattices are used [22-24]. Ref. [23] quotes results from two lattice 
spacings. In Table [2] we present the results from their finer lattice. To determine the strange 
quark mass parameter k 8 , we interpolate the vector meson mass to the physical <fi meson mass, 
= 1.01946(19) GeV. We find k s = 0.13502(6). This is our "method I" for determining 
the k value corresponding to the strange quark mass. Using Eq. (j4]) and setting m K to the 
experimental value for (m 2 K+ +m 2 Ko )/2 gives a value in very close agreement: k s = 0.134981(9). 

The raw data for the heavy-light meson masses are given in Table |3j To find the physical 
values of the heavy-light meson masses, we extrapolate for each heavy quark mass linearly in 
m q , see Eq. (JSJ). The fits are shown in Fig. [TJ The quark mass dependence is linear to very 
good accuracy. This is in contrast to the findings of, e.g., Ref. [25]. In the final step, the 
calculation of the decay constants, the physical values of the c and b quark masses will be 
reached by interpolating or extrapolating the heavy-light meson mass to the D or B mass and 
the heavy-strange meson mass to the D s or B s mass. 

In a quenched calculation, different methods to choose the input for determining physical 
parameters may give different answers. In order to investigate the influence of this arbitrariness 
we also use the heavy-light spectrum to determine k s and call this procedure "method II" . We 
consider the splitting between mesons with a heavy quark and a strange quark (generically 
denoted by M s ) and a meson with a heavy quark and a quark with the u, d quark mass (denoted 
by Mi). In our data, as well as in experiment, the M s — Mi mass difference is fairly independent 
of the heavy quark mass. To fix k s in method II, we choose a heavy quark close to the charm 
mass from our simulation points, namely k = 0.129, and set the splitting between the M s and 
the Mi masses equal to the experimental value for the D meson, vn^ s — mp = 98.85(30) MeV. 
The corresponding value for the strange hopping parameter is k s = 0.134929(15). 

We calculate the pseudoscalar decay constants from the improved axial vector current A 1 ^ 

A{ = Z A (1 + ab A m q ) (A 4 + c A ad 4 P) , (6) 

where A^(x) = q\ x l^q2x and P(x) = q^b^x- We take the nonperturbatively determined 
values for Z A from [26] and for c A from [18]. For our calculation, this gives Z A = 0.8338 and 
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Ref. 


a 1 [GeVJ 


quark action 


gauge action 


m p [Ge V] 


this work 


4.97 


clover 


Wilson 


0.849(38) 


[21] 


cont 


clover 


Wilson 


0.797(13) 


[22] 


1.33 


chirally improved 


Liischer-Weisz 


0.791(42) 


[22] 


1.29 


fixed point 


fixed point 


0.828(25) 


[23] 


2.09 


overlap 


Liischer-Weisz 


0.79(2) 


[24] 


1.60 


overlap hypercube 


Wilson 


1.017(40) 



Table 2: p meson masses from quenched lattice calculations. The lattice scale has been determined 
using tq = 0.5 fm in all calculations except in [23] where ro = 0.56 fm is used. The quoted errors are 
only statistical. 

ca = —0.01967. The coefficient bA is calculated from 1-loop perturbation theory [27]. Using a 
boosted coupling g% — > g^/uQ with u = (^TrUp) 1 ^ 4 , we find &4 = 1.2143 which is close to the 
result one finds using the tadpole-improved scheme of [28]. A non-perturbative determination 
of bA on coarser lattices (J3 < 6.4) [29] also gives values in agreement with boosted perturbation 
theory within errors. 

The meson matrix elements of the currents 

/ (0) = ^(0|A|M), 

/W = ^(°l a ^P|M) = -^sinh(aM)(0|P|M), 

/ = ±(0\A{\M), (7) 
are related to the amplitudes by 

(0) = -O.f^Z ^PAA (ox 



y/MVA% S p 



f» - 2^sinh( M)-^, (9) 

where M denotes the meson mass. The convention for the factors of a/2 corresponds to the 
normalization where f n ~ 130 MeV. 

3. The fit of the light meson decay constants according to Eq. (jHJ) is shown in Fig. [2j Mesons 
with degenerate and nondegenerate quark masses fall on the same straight line. For f n , the 
value at the physical u, d quark mass, and fx, the value extrapolated to the averaged strange 
and u, d quark mass, we find 

U = 140(4) MeV, f K = 153(4) MeV. (10) 

This result for f n agrees well with the value of 137(2) MeV determined by [21] using an ex- 
trapolation to the continuum from coarser lattices. Both values are slightly larger than the 
experimental value of f w + = 130.7(4) MeV. Our value for fx is 6% or 2a lower than the result 
from [21] of 163(1) MeV. The experimental value is f K + = 159.8(15) MeV. 
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amps 


af° 


af 1 


af 


0.11500 


0.13519 


0.8363(15) 


0.0371(11) 


0.0532(17) 


0.0423(13) 


0.12100 


0.13519 


0.6676(13) 


0.0417(14) 


0.0496(17) 


0.0432(14) 


0.12900 


0.13519 


0.4065(11) 


0.0475(13) 


0.0417(13) 


0.0435(12) 


0.13000 


0.13519 


0.3685(12) 


0.0478(13) 


0.0399(13) 


0.0431(12) 


0.11500 


0.13498 


0.8431(12) 


0.0383(12) 


0.0551(19) 


0.0437(13) 


0.12100 


0.13498 


0.6747(11) 


0.0429(12) 


0.0517(16) 


0.0446(12) 


0.12900 


0.13498 


0.4145(10) 


0.0488(15) 


0.0431(14) 


0.0448(14) 


0.13000 


0.13498 


0.3765(09) 


0.0490(13) 


0.0412(12) 


0.0443(12) 


0.11500 


0.13472 


0.8517(11) 


0.0402(12) 


0.0584(19) 


0.0460(13) 


0.12100 


0.13472 


0.6836(10) 


0.0446(13) 


0.0541(16) 


0.0466(14) 


0.12900 


0.13472 


0.4242(08) 


0.0508(13) 


0.0453(13) 


0.0469(12) 


0.13000 


0.13472 


0.3866(08) 


0.0507(13) 


0.0430(12) 


0.0460(12) 



Table 3: Pseudoscalar heavy-light meson masses and decay constants at the simulation points. 




Figure 1: Chiral extrapolation of meson masses. On the left, light vector meson masses, on the 
right, heavy-light pseudoscalar meson masses for the heavy hopping parameters k = 0.115 (circles), 
k = 0.121 (squares), k = 0.129 (diamonds) and k = 0.130 (triangles). Open symbols denote the 
simulation points, closed symbols the chiral extrapolation. 

The SU(3) flavor breaking ratio of the light decay constants in our calculation turns out to 
be relatively small. We find 

f K /U - 1 = 0.088(12) . (11) 

Our number is substantially lower than the experimental value of 0.222, but is consistent with a 
recent quenched calculation using overlap fermions [30], which finds /k/Jtt — 1 — 0.09(4) using 
the same scale setting with r = 0.5 fm. It is also consistent with other quenched determinations 
(see [31]). 

4. Next we consider the heavy-light decay constants. To determine values at the physical 
quark masses, we extrapolate or interpolate the decay constants separately in the light and the 
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am 

q 

Figure 2: Chiral fit of fight meson decay constants. The chiralfy extrapofated vafue is denoted by the 
filled circie. 

heavy quark mass. For the fits in the light quark mass we use a function of the form (jSj) with 
m q being the mass of the quark with the light k parameter of the heavy-light meson instead of 
the average quark mass. For the extrapolation to the b quark mass and also for an interpolation 
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Figure 3: On the left, chiral extrapolation of heavy-light decay matrix elements. Symbols have the 
same meaning as in the right part of Fig.[TJ On the right, heavy quark mass dependence of heavydight 
decay matrix elements. Squares denote strange, and diamonds denote physical light quarks. Closed 
symbols denote heavy quark masses extrapolated to the b or interpolated to the c quark mass. 

to the c quark mass we use a formula motivated by HQET (see e.g. [32]). In the heavy quark 
limit, matching of the decay matrix element in the effective theory to the matrix element in full 
QCD introduces logarithmic corrections in the heavy quark mass which have to be resummed. 
In addition, power corrections in 1/ tuq have to be added. Since the extrapolation to the b mass 
in our case is rather short, the precise form of the extrapolation formula is not important. We 
use only the lowest order running for a s , and take the heavy-light meson mass (and not the 
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quark mass tuq) as an expansion (scale) parameter: 



Here 70 = —4 is the leading order anomalous dimension of the axial vector current, and bo — 11 
is the leading coefficient of the QCD (3 function for zero dynamical flavors. The fits and the 
interpolated values are shown in Fig. [31 The values of the fit parameters are c = 0.55(4) 
GeV 3 / 2 , ci = -0.66(19) GeV and c 2 = 0.38(21) GeV 2 if the light quark mass is the u, d quark 
mass, and c = 0.59(4) GeV 3 / 2 , c x = -0.70(14) GeV and c 2 = 0.39(15) GeV 2 for the s quark. 



Decay constant ratios 


fD a l 'Id 


IbJ Jb 


fDs 1 /b 3 


Id/ /b 


1.068(18)(20) 


1.080(28)(31) 


1.069(28)(160) 


1.082(42)(168) 



Table 4: Ratios of heavy-light decay constants. The first error is statistical, and the second systematic. 
The systematic errors are discussed in the text. 

Our final results for the ratios of heavy-light decay constants are presented in Table HJ and 
the heavy-light decay constants are given along with a comparison in Table [5j 

Estimation of systematic errors is notoriously difficult. One source of uncertainty concerns 
setting the quark masses to their physical values. For the strange quark this can be estimated 
by comparing the results from our methods I and II and suggests an error of 4 MeV for fu s 
and fs s - For the u, d quarks a chiral extrapolation is required. The corresponding error is 
difficult to estimate. Our data are consistent with the simplest linear chiral extrapolation. 
Quenched chiral perturbation theory provides a more sophisticated formula. However, it is not 
clear if it is applicable to our data. The uncertainty in fixing the heavy quark mass can be 
estimated by comparing the difference between the mass fixed from quarkonium and from the 
heavy-light meson system. Since the r\ c meson mass using the charm quark hopping parameter 
determined from the D s meson agrees with the physical value, we assume that this uncertainty 
is rather small in our calculation. In addition, for the B system there is an uncertainty from the 
extrapolation in the heavy quark mass. The difference between a quadratic fit to the matrix 
elements fVM and a quadratic fit to $ is very small and changes the values for the decay 
constants by less than 1 MeV. If only the three lighter heavy quark masses are included in the 
extrapolation to the b mass, fs a (/b) changes by —3 (+1) MeV. 

Since we have results only from one lattice spacing, we cannot perform a continuum ex- 
trapolation from our data alone and have to estimate the discretization effects as a systematic 
error. Leading discretization effects are 0(a 2 ). A rough estimate of them can be obtained by 
squaring the 0(a) corrections appearing in the Symanzik improvement program. For the charm 
quark, the correction proportional to ca is small, around 2%, while the term proportional to 
the quark mass and bA is around 10% of the size of the matrix element itself. The square of 
the sum of these variations is around 1%, which we take as our estimate for the discretization 
error of fo and fo s ■ A similar consideration for the B and B s systems results in an estimate of 
a discretization error of roughly 12%. For the error in the renormalization constants we use the 
estimate given for Za in Ref. [26] of 1%. Since the heavy-light meson masses in lattice units 



7 



0.6 




2 °- ; 



6 



6 




Hi 



0.3 



0.4 



0.5 



l/M[GeV] 



Figure 4: Difference of the kinetic mass and the rest mass for heavy-light mesons with the light 
hopping parameter n = 0.13498. The line at zero is plotted to guide the eye. 

in our simulation increase up to values of ~ 0.8 one might be concerned about cutoff effects in 
the dispersion relation for the heavy-light meson. We therefore compare the kinetic mass [7] 
Mkm calculated from 



where M is the rest energy of the meson. Results for the light quark mass close to the strange 
quark mass are shown in Fig. HI We find that discretization errors in the dispersion relation 
are smaller than the statistical errors. 

The finite volume effects of the ratio Jbs/Jb have been investigated in the framework of 
heavy meson chiral perturbation theory [33]. For quenched lattices of spatial extent 1.6 fm and 
pseudoscalar meson masses around 500 MeV (which corresponds to the smallest quark mass 
used in our simulation) they are quoted to be around 1% or smaller. It is plausible that the 
finite volume effects for D and D s mesons are of similar size. 

We estimate the total systematic error due to discretization effects, errors in Z4, finite volume 
effects and ambiguities in fixing the physical quark masses by collecting all contributions and 
adding them in quadrature. It is given as the second error in Table [51 

The ratios are less sensitive to some of the systematic effects. The dominating ones for 
/d s //b 8 and /d//b are the discretization effects. For /d s / Id and fnj fB we find a variation 
depending on how the strange quark mass is set, while the estimated discretization effects are 
smaller than the statistical errors. The uncertainties from fixing the physical quark masses and 
the discretization errors (added in quadrature) are given as the second error in Table HI 

We have used the value of the Sommer parameter r$ = 0.5 fm to set the scale in physical 
units. This choice allows for a direct comparison with previous lattice determinations (see 
below) but is not universally accepted. With a different value of the Sommer parameter our 
results have to be modified accordingly. The variation of the decay constants if r is changed 
by ±10% is given as third error bar for our results in Table 00. 

5. Finally, we compare our results to other lattice calculations of decay constants. There exist 
recent quenched results for fu a from nonperturbatively 0(a) improved clover fermions [34-36] 

^"We note, however, that a value of r Q = 0.45 fm leads to seemingly unphysical results. In particular the 
SU(3) breaking in the meson masses and decay constants becomes very small. Also, different methods to set 
the strange quark mass produce more noticeable differences in the results. 
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Ref. 


Nf, HQ action, scale 


/rUMeV] 


/d [MeV] 


Lattice 


this work 
[34] 
[35] 
[36] 

[9] 

[25] 

[42] 


0, clover, ro = 0.5 fm 

0, clover, ro = 0.5 fm 

0, clover, r = 0.5 fm 

0, clover, ro = 0.5 fm 

0, mod. Fermilab, ro = 0.5 fm 

0, overlap, f n 

2 + 1, Fermilab, T spectrum 


220(6)(5)(11) 

243(2) (§3) 

252(9) 

225(6) 

237(5) 

266(10)(18) 

249(3) (16) 


206(6)(3)(22) 
222(3)(°|) 

235(8) (14) 
201(3)(17) 


Experiment 


[3] 
[4] 
[2] 




280(12)(6) 
283(17)(16) 


223(17)(3) 


Ref. 


Nf, HQ action, scale 


/iUMeV] 


/b [MeV] 


Lattice 


this work 

[34] 

[38] 

[39] 
[12] 
[13] 


0, clover, ro = 0.5 fm 

0, clover, r = 0.5 fm 

0, clover+static, r = 0.5 fm 

0, clover+static, r = 0.5 fm 

2 + 1, NRQCD, T spectrum 

2 + 1, NRQCD, T spectrum 


205(7)(26)(17) 
240(4)(2) 
205(12) 
191(6) 
260(7) (28) 


190(8)(23)(25) 
217(5)(Jg) 

216(9)(20) 


Experiment 


[1] 
[41] 


experiment 
UTfit 


227(9) 


229(g) (|?) 



Table 5: Heavy- light decay constants from lattice calculations and experiment for the D (upper table) 
and for the B system (lower table). For the lattice calculations, the number of flavors in the simulation 
(Nf), the heavy quark (HQ) action, and the quantity used to set the scale are also indicated. The first 
error bar is the statistical, the second (where given) the systematic error except for the uncertainty in 
ro- For our work we quote a third error assuming a ±10% uncertainty in the physical value of ro. For 
the result from [36] we quote the value from the finest lattice instead of the continuum extrapolated 
result. 

for a range of lattice spacings (0.03 < a < 0.1 fm) as well as for overlap quarks [25]. The 
comparison with the clover results is particularly interesting because it sheds some light on the 
discretization effects and might indicate the possibility of a joint continuum extrapolation. We 
plot the clover data in Fig. [5] as a function of the squared lattice spacing together with the 
overlap data. First we notice that on coarser lattices there is a discrepancy between the clover 
data of Refs. [34] and [35]. The discrepancy corresponds roughly to the difference one obtains 
when ca values from different nonperturbative calculations for a meson mass > 2.4 GeV are 
used, as discussed in [34]. Furthermore, the work [35] uses a nonperturbatively determined 
value for b& [37]. On the finer lattice of Ref. [34] ((3 = 6.2) the value used in Ref. [35] is about 
6% larger than the perturbative number, which according to our estimates would affect the 
decay constants by at most 2%. At f3 = 6.0 the difference is even smaller. On the finest lattice 
used by [35] (f3 = 6.45) the difference between the perturbative and nonperturbative values of 
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\>a is ~ 7%, which translates on a fine lattice into only a very small difference in the decay 
constants. In a more recent calculation [29] the nonperturbative value at that (3 value has come 
into agreement with perturbation theory, as mentioned in Section 2. 

Our data is in good agreement with the value obtained by Jiittner on his finest lattice [36]. 
The overlap value from [25] is on the other hand substantially larger. Being determined on a 
relatively coarse lattice it might be affected by discretization errors. It is important that all 
data shown in Fig. [5] come from lattices with similar spatial extent between 1.5 and 1.6 fm. So, 
finite size effects can be expected to be roughly the same in all calculations. 



280 r , 1 r 




,1 1 1 1 1 1 1 1 1 1 1 

0,05 0,1 0,15 0,2 0,25 

a 2 [GeV" 2 ] 



Figure 5: Lattice spacing dependence of quenched fn a from O(a) improved clover quarks (this work, 
star), (UKQCD [34], diamonds), (ALPHA [35], squares), (Jiittner [36], circle), and overlap quarks 
(Ref. [25], triangle). The error bars show statistical and fitting uncertainties only. The scale is set 
using ro = 0.5 fm with the exception of [25] where f n is employed for the conversion of the decay 
constant to physical units. If ro = 0.5 fm is used instead, their lattice spacing decreases by 12%, which 
would increase their result for the decay constant even further. 

Our result for fs s is consistent with the quenched calculations of [38,39], but considerably 
lower than the nonrelativistic (but unquenched) calculation of [12]. The fit to the standard 
model gives a value with a relatively small error in between these two numbers. 

For fjj and fs, the values obtained from lattice calculations are consistent with the experi- 
mental results. Since the experimental errors are still large, this comparison is not conclusive, 
however. 

6. Let us summarize our main findings. We have calculated decay constants of heavydight 
pseudoscalar mesons on a very fine quenched lattice using clover fermions. Our extrapolations 
to the b quark mass appear reasonable. Nevertheless, from a comparison of the results at the 
charm mass to data obtained on coarser lattices we obtain the impression that discretization 
errors with the relativistic formalism adopted here are still significant for the b sector, unless 
the inverse lattice spacing becomes larger than ~ 10 GeV. 

Our results and those of Ref. [36] for f Da are 10 — 15% smaller than the central values quoted 
for other recent lattice calculations, and roughly 20% smaller than recent experimental values. 
Eventually one would like to determine the decay constants to an accuracy of a few percent. 
Our work and the result of Ref. [36] indicate that discretization errors for the clover results on 
lattices with a" 1 < 2 — 3 GeV are too large to reach this precision, and that even a continuum 
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extrapolation from a set of coarser lattices has a large uncertainty for heavy quarks. On the 
other hand, we do not find any source of large systematic errors, other than quenching, that 
could affect our calculation. It seems, therefore, that the new lattice results on fine lattices (this 
work and [36]) indicate a relatively small value for fu s from lattice QCD. Quenching effects are 
notoriously difficult to estimate. However, since in previous calculations with a" 1 ~ 2 GeV [15] 
it was found that the quenching error is insignificant with our choice of lattice parameters, we 
expect that they will not be too large. 

The systematic uncertainties on our results are larger for the B system than for the D 
system and more difficult to estimate reliably. Our results are in agreement with several other 
recent lattice calculations, but smaller than the values from recent unquenched calculations 
using nonrelativistic methods. 

We find a rather small 577(3) symmetry breaking ratio of the heavy-light and light decay 
constants compared to experiment and also to several recent unquenched lattice calculations. 
The difference between our numbers and the unquenched results may be partially due to the 
use (see e.g. [42]) of a chiral extrapolation formula for the unquenched data which is inspired 
by chiral perturbation theory and predicts a particularly strong decrease of the decay constant 
at lighter quark mass values than are accessible in the simulation. This is in contrast to the 
use of a simple linear extrapolation in our calculation. 
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